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(Dated: November 12, 2010) 

A search for sidereal modulation in the flux of atmospheric muon neutrinos in IceCube was 
performed. Such a signal could be an indication of Lorentz-violating physics. Neutrino oscillation 
models, derivable from extensions to the Standard Model, allow for neutrino oscillations that depend 
on the neutrino's direction of propagation. No such direction-dependent variation was found. A 
discrete Fourier transform method was used to constrain the Lorentz and CPT- violating coefficients 
in one of these models. Due to the unique high energy reach of IceCube, it was possible to improve 
constraints on certain Lorentz-violating oscillations by three orders of magnitude with respect to 
limits set by other experiments. 

PACS numbers: 11.30.Cp,04.60.-m,14.60.St,95.55.Vj 



I. INTRODUCTION 

Lorentz invariance and CPT symmetry, which com- 
bines charge conjugation (C), coordinate reflection (P), 
and time reversal (T), are fundamental symmetries of 
quantum field theory (QFT). To date, no experimental 
evidence for a violation of either symmetry has been ob- 
served, despite a wide variety of experimental investiga- 
tions 1 3]. However, it remains worthwhile to continue 
to test these fundamental symmetries, with different ex- 
periments and different types of particles, at higher en- 
ergy scales or with improved sensitivity. Observation of a 
violation of one of these symmetries would be an indica- 
tion of new physics, and possibly point the way towards 
a unified theory or a theory of quantum gravity. 

There is reasonable motivation to expect that Lorentz 
invariance and CPT symmetry do not hold all the way 
to the Planck scale (Mp « 10 19 GeV), due to a discrete 
structure of spacetime or interactions with a spacetime 
foam [4, 5], for example. Neutrinos are sensitive probes of 
possible low energy effects of the breaking of these sym- 



metries, because they have very high Lorentz factors and 
they do not interact by the strong or electromagnetic 
forces. Signatures of Lorentz and CPT-violating pro- 
cesses in the neutrino sector may include oscillations with 
unique energy dependencies, direction-dependent oscil- 
lations that violate rotational invariance, or deviations 
from the anticipated behavior based on the L/E ratio of 
the experiment @. 

The Standard Model (SM) of particle physics is be- 
lieved to be the low-energy limit of a more fundamental 
theory. Such an extension of the SM is typically assumed 
to unite QFT and General Relativity at the Planck scale, 
and to provide a coherent theory of quantum gravity. 
To look for signatures of quantum gravity without this 
fundamental theory, a phenomenological description is 
necessary. The Standard Model Extension (SME) [7] is 
an effective-field-theory framework that provides such a 
phenomenological description at experimentally accessi- 
ble energies, and has guided numerous searches for pos- 
sible signatures of Lorentz invariance violation and CPT 
violation 0]. 
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The IceCube detector @, located at the South Pole, is 
designed for detecting astrophysical neutrinos of all three 
flavors. Due to its unique size, it has an unprecedented 
event rate for high energy atmospheric neutrinos. Data 
taken while IceCube operated in a partially-completed, 
40-string configuration, was used to search for a peri- 
odic variation as a function of right ascension, a possible 
consequence of a Lorentz-violating preferred frame. A 
discrete Fourier transform (DFT) method was used to 
constrain Lorentz and CPT violating coefficients in the 
SME, in the context of a direction-dependent neutrino 
oscillation model that violates rotational invariance. 



II. THE VECTOR MODEL 

The SME adds to the SM Lagrangian all terms that 
can be constructed with SM and gravitational fields, but 
that may also violate Lorentz or CPT symmetries. The 
coefficients for these processes have Lorentz indices and 
represent background tensor fields. Physically observable 
phenomena depend on contractions between these tenso- 
rial coefficients and the particle momentum. A subset of 
the SME, known as the "minimal" SME 6], includes all 
observer- independent, renormalizable, Lorentz and CPT- 
violating processes. Energy and momentum are still con- 
served, and spin-statistics and gauge invariance are pre- 
served. Right-handed neutrinos are still assumed to de- 
couple and remain undetectable. Neutrino masses are 
treated the same as in the SM. 

The effective Hamiltonian from the minimal SME, de- 
scribing Lorentz-violating oscillations between neutrino 
flavor states a and b, is 16!] 

OcffU = [(a L )abPn ~ ( c L)ab Pf*P»] > (!) 

where E is the neutrino energy, the neutrino four- 
momentum, and /i and v are Lorentz indices. The coeffi- 
cients (a-L)ab nave mass dimension 1 and lead to Lorentz- 
violating and CPT-violating interactions. The coeffi- 
cients {cl) 1 ^ nave mass dimension and lead to Lorentz- 
violating interactions. After some approximations appli- 
cable to the length and energy scale of atmospheric neu- 
trinos, as well as some assumptions about which com- 
ponents of the interaction tensors are non-zero, a subset 
model known as the "vector model" can be derived [6j. 
The vector model is convenient for studying possible side- 
real variations in the atmospheric neutrino flux. 

In the vector model, the only non-zero com- 
ponents of the interaction tensors in Eqn. [T] are 

( a £V ) ( a L) M T ' ( C i)^r > and ( c £) M T > a11 assumed to be 
real. Only vacuum oscillations between neutrino flavor 
states /i and r are included. These assumptions are made 
in a Sun-centered celestial-equatorial coordinate system. 
The z-axis is aligned with the Earth's rotational axis 
and the x-axis points towards the vernal equinox. Mass- 
induced oscillations between v„ and v r are not included 
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FIG. 1. Simulation of the sinusoidal signal predicted by 
Eqn. [2] with a x = 2 x 1CT 23 GeV. 



in the vector model. However, they were included in the 
simulation of the expected neutrino flux. 
The Ufx survival probability is then Q 



in 2 {L {A s ) fiT sm(a + tp ) 
+ ( A c) MT cos {a + <p ) ), 



(2) 



where L is the propagation distance, a is the neutrino's 
right ascension, and ifQ is an arbitrary phase offset. Drop- 
ping the flavor subscripts, A s and A c are defined as 



A. = N Y (af - 2Ecl x ) - N* (a Y L - 2Ec TY ) 



A r 



-N x (ai 



2EcD 



N Y ( 



2Ecl Y ) 



(3) 
(4) 



The survival probability for antincutrinos, P p _>£> , is 
given by changing the sign of the az, coefficients. The 
oscillation probability depends intrinsically on the direc- 
tion that the neutrino propagates through space, violat- 
ing rotational invariance. The N X ( Y *> are unit propaga- 
tion vectors for the neutrino: 



N x = sin (9) cos (<p) 
N Y = sin (9) sin (<p) , 



(5) 
(6) 



where 9 — ir/2 + 5, if = n + a, and S is the declination 
of the incident neutrino. Figure [1] shows an example of 
the anticipated sinusoidal signal in IceCube, as predicted 
by Eqn. [1 with a x = 2 x I CT 23 GeV and the detector 
configuration and live time discussed in the next section. 



III. THE EVENT SAMPLE 

When completed in 2011, IceCube [8j, [9( will consist of 
86 strings. Each string includes 60 digital optical mod- 
ules (DOM), for a total of 5,160 DOMs. A DOM is a 
single photomultiplicr tube and associated electronics in 
a glass pressure sphere [Io| . The instrumented part of the 
array extends from 1450 m to 2450 m below the surface 
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of the ice. Horizontally, 78 of the strings are 125 m apart 
and spread out in a triangular pattern over a square kilo- 
meter. Vertical DOM spacing is a uniform 17 m for these 
78 strings. A subset of the detector, known as "Deep- 
Core" , consists of 8 specialized and closely spaced strings 
of sensors located around the center IceCube string. 

This analysis used data from 359 days of live time while 
operating in a 40-string configuration, from April 2008, 
to May 2009. No DeepCore strings had been installed at 
that time. The event sample is a subset of the data used 
for an unfolding of the atmospheric muon neutrino spec- 
trum [llj ]. Triggering, filtering, and background rejection 
are discussed in detail in [Tl[ . The 40-string detector was 
roughly twice as long in one horizontal direction as in the 
other. However, this azimuthal dependence of the detec- 
tor shape conveniently canceled out due to the sidereal 
rotation of the Earth (and thus, the detector around its 
vertical axis). 

IceCube detects the Cherenkov radiation from charged 
particles produced in charged current (CC) and neutral 
current (NC) interactions between incident neutrinos and 
nucleons in the ice. If the incident neutrino is a or v^, 
a muon or antimuon is produced and undergoes radia- 
tive energy losses as it propagates, creating additional 
Cherenkov radiation. The muon directions are recon- 
structed from records of photon arrival times at DOMs 
participating in each event. The mean angular devia- 
tion between the direction of the parent neutrino and 
the muon is less than a degree for the energy range of 
this analysis. Additionally, simulation and reconstruc- 
tion studies indicate that muon angular resolution is typ- 
ically between 0.5° and 1°, depending on the angle of in- 
cidence and the muon energy. Hence, the reconstructed 
muon direction provides a good estimate for the neutrino 
direction. 

Background events in the data were down-going at- 
mospheric muons, or coincident muons, that were recon- 
structed as up-going events. Rejection of this background 
was done in several stages, beginning with triggering and 
local coincidence checks on the DOMs 13], and software- 
based filtering at the South Pole [llj . Then, before 
more computationally intensive reconstructions were per- 
formed during off-line processing, unusable events were 
removed by selection cuts based on zenith angle and track 
quality parameters. Finally, using boosted decision trees 
(BDT) [LJ] , we obtained a sample of 7882 muon neutrino 
events in the zenith range 97° to 120°, with negligible 
background. Background contamination was estimated 
to be less than 1%, based on testing the BDTs with sim- 
ulated atmospheric muon and neutrino data sets. This 
value was then verified by comparing the data passing 
rate as a function of BDT cut value to the predicted rate 
from atmospheric muon and neutrino simulation. These 
event selection cuts also eliminated localized events from 
electromagnetic showers induced by v e CC interactions 
and hadronic showers due to NC interactions. 

As discussed in [llj] , there was a statistically significant 
excess of events in data (or deficit in simulation) in the 
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FIG. 2. RA distribution of events in data. Vertical error 
bars are statistical uncertainty only. Fluctuations in the data, 
above and below the mean (horizontal line), are consistent 
with statistical variations, \ 2 P er bi n f° r a straight-line fit to 
the mean is 0.9. 

zenith region 90° to 97°, the origin of which could not be 
verified at the time of the analysis. Hence, that region 
was not used. The vector model we adopted assumes 
that only real components belonging to the plane per- 
pendicular to the Earth's axis are non-zero and ignores 
any coupling between the z-component of the neutrino 
momentum and the Lorentz violating coefficients of the 
SME Q. By considering only events in the zenith re- 
gion from 97° to 120°, where the x and y components of 
the neutrino momenta dominate, the impact of this ar- 
bitrary assumption on how the Lorentz violating field is 
aligned with respect to our preferred coordinate system 
is minimized. 

Figure [2] shows the distribution in right ascension (RA) 
for events in the data. This histogram has 32 bins from 
to 360° in RA, the same binning that was used to com- 
pute power spectral densities with DFTs as discussed in 
the next section. We estimate from simulation that about 
90% of the events are from atmospheric neutrinos in the 
energy range 200 GeV to 13 TeV, and 99% from atmo- 
spheric neutrinos in the energy range 100 GeV to 55 TeV. 



IV. METHODOLOGY AND RESULTS 

The DFT analysis methodology was adapted from 
Ref. [TH, where the MINOS collaboration looked for side- 
real variations in the NuMI beam line, using the MINOS 
near detector. DFTs were computed using FFTW [l6[ in 
the ROOT [l7j framework. Under the vector model, the 
muon neutrino survival probability varies with RA with a 
modulation frequency of 4w0, where = 27r/23h56 min 
is the Earth's sidereal frequency. So we are interested in 
the n = 4 mode of a DFT of the event rate as a function 
of RA. 

First, the data were checked for consistency with the 
hypothesis of no sidereal signal. For each of 10,000 trials, 
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the RA of all data events were randomized and the power 
spectral densities (PSD) in modes n = 1 through n = 4 
of a DFT were computed. The PSDs for the true data 
were then computed and compared to these "noise-only" 
distributions. The data was consistent with no signal in 
any of the modes. In particular, the PSDs for data were 
less than 34% of the noise-only trials for n = 1, 92% for 
n = 2, 31% for n = 3, and 98% for n = 4. 

With the absence of a sidereal signal, we were able 
to set upper limits on the SME coefficients in the vec- 
tor model. The flux models of [la , Il9f were assumed 
for conventional and prompt atmospheric neutrinos. The 
predictions for neutrinos from pions and kaons were ex- 
tended to higher energies by fitting a physics-motivated 
analytical equation based on energy and zenith angle 
( [20] | and chapter 7 of [ll|), in an overlapping region with 
the detailed calculations of [Hj . 

Normally, we could ignore ^-induced muons for an 
atmospheric neutrino analysis. However, if some of the 
Vfj, and are oscillating to v T and v T according to the 
model that we want to test, then that flux of v T and 
v T has to be accounted for. A v T (D T ) could undergo 
a CC interaction in or near the detector, producing a 
tau lcpton, which could then decay into (among other 
things) a muon (branching ratio about 17%). Detection 
of these muons would dampen the signal we are looking 
for, i.e. a disappearance of muon neutrinos. This effect 
was accounted for through toy Monte Carlo studies and 
v r (v T ) simulation. About 6% of the events lost due to 
oscillations induced by the coefficients, and about 9% 
of the events lost due to oscillations induced by the cl 
coefficients, are recovered. The difference between the 
two cases is due to the fact that the mean energy of 
affected events is higher in the case of the Cl coefficients, 
and detector efficiency increases with energy. 

In each of 400 toy Monte Carlo (MC) experiments, sim- 
ulated events were drawn from a distribution matching 
the energy and zenith dependence of atmospheric neu- 
trinos. RA's were randomly assigned to the simulated 
events in each trial. The physics parameters of the vector 
model were incrementally increased, and the simulated 
events re-weighted according to their survival probabil- 
ity under the vector model, until a PSD greater than 
the 99.87 percentile (equivalent to a 3-sigma threshold) 
of the PSDs from the 10,000 noise-only toy experiments 
was obtained. The values found in each of these 400 
trials were then averaged to estimate the sensitivity to 
a sidereal signal described by the vector model. These 
trials were performed independently for each coefficient. 
While adjusting one coefficient, the other three were held 
at zero. 

The sensitivity depends on the zenith and energy dis- 
tribution of atmospheric neutrinos, and is thus affected 
by uncertainties in these quantities. Theoretical and ex- 
perimental uncertainties in the zenith distribution are 
small, a 3% uncertainty in the predicted ratio of the ver- 
tical to horizontal atmospheric neutrino flux 11811 . and 
angular resolution on the order of a degree [llj . Toy 



MC experiments with the simulated zenith distribution 
modified according to these uncertainties showed that the 
impact on the sensitivity is negligible. 

Uncertainties in the energy distribution do not affect 

X( Y) 

sensitivity to oscillations driven by the a L coefficients. 
However, uncertainties in the spectral index for atmo- 
spheric neutrinos [l8l [l9j . and uncertainties in the en- 
ergy dependence of the detector efficiency, both affect 
sensitivity to the c^ x ^ TY ^ coefficients. Uncertainty in 
the spectral index is primarily due to uncertainty in the 
energy distribution of the cosmic ray flux [12, . The 
uncertainty in the spectral slope of the proton component 
of the cosmic ray flux is assumed to be ±0.03, and for 
the helium component (which makes up roughly 30% of 
the cosmic rays in the energy region of this analysis) it is 
assumed to be ±0.07. Combining these two factors, after 
scaling them by the fraction of their contribution to the 
total flux, leads to an estimated ±0.05 uncertainty in the 
spectral index. Toy MC experiments with simulated at- 
mospheric neutrino events re- weighted to account for this 
uncertainty in the spectral index showed a ±7% change 
to the sensitivity for the c£ T coefficients. 

DOM sensitivity and ice property uncertainties affect 
the energy dependence of the detector's effective area, 
and hence the distribution of neutrino energies for events 
in the data. Specialized simulated data sets with ±10% 
enhanced photon populations were used to estimate this 
uncertainty, in a manner similar to the evaluation of 
DOM sensitivity and ice property uncertainties for the 
unfolding analysis discussed in [llj . A ±10% change in 
the number of photons observed by the DOMs leads to 
a ±0.05 change in the spectral index for the energy dis- 
tribution of detected neutrinos, which in turn leads to a 
±7% change in sensitivity for the ( ^ X ^ TY ^ coefficients. 

The uncertainty in DOM sensitivity was estimated to 
be ±8%, based on the measured uncertainty in PMT 
sensitivity [l(i| |. This was directly scaled to the ±10% 
change in the number of simulated photons striking the 
DOMs. The average photon flux was estimated to change 
by ±12%, due to uncertainties in scattering and absorp- 
tion, using a diffuse flux approximation [llj . Added in 
quadrature, DOM sensitivity and ice property uncertain- 
ties lead to a ±15% uncertainty in the number of detected 
photons, and a ±11% uncertainty in sensitivity to the 



TX(TY) 



coefficients. 



The following upper limits have been set on the SME 
coefficients, at the 3 sigma level: 



,X „Y 



'L > 



a' L < 1.8 x 10 



-23 



GeV, 



and 



c[ x ,cl Y < 3.7 x 10 



-27 



(7) 



(8) 



A net 13% increase (0.4 x 10 27 ) has been added to the 
upper limit for c^ x ^ TY \ to account for systematic flux 
(7%) and detector (11%) uncertainties added in quadra- 
ture. 
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V. CONCLUSION 

We have found no sidereal variation in the atmospheric 
muon neutrino event rate in IceCube. In the context of 
the SME, we found no evidence for a violation of Lorentz 
or CPT symmetries due to a preferred reference frame. 
The LSND [H and MINOS Qj| HH Collaborations also 
did not see sidereal variations in the number of neutrinos 
detected from their respective beam lines. With their far 
detector [H, MINOS found af < 5.9 x 1CT 23 GeV and 
a\ < 6.1 x 1CT 23 GeV, c\ x and c£ Y < 0.5 x 1CT 23 . Our 
results of a x ,a\ < 1.8 x 1CT 23 GeV, and c T L x ,cl Y < 
3.7 x 1CV 27 , have improved upon these limits by a factor 



of three for the 



coefficients and by three orders of 



magnitude for the c L y ' coefficients, due to the long 
baseline of atmospheric neutrinos and the high energy 
reach of IceCube. 
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